Entanglement at the boundary of spin chains near a quantum critical point and in 

systems with boundary critical points 
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We analyze the entanglement properties of spins (qubits) attached to the boundary of spin chains 
near quantum critical points, or to dissipative environments, near a boundary critical point, such 
as Kondo-like systems or the dissipative two level system. In the first case, we show that the 
properties of the entanglement are significantly different from those for bulk spins. The infiuence of 
the proximity to a transition is less marked at the boundary. In the second case, our results indicate 
that the entanglement changes abruptly at the point where coherent quantum oscillations cease to 
exist. The phase transition modifies significantly less the entanglement. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Hk 



I. INTRODUCTION 

Quantum phase transitions have attracted intense re- 
search activities on various fields of physics^. Whereas 
classical phase transitions are driven by thermal fluctu- 
ations, quantum transitions are induced by a parameter 
which enhances quantum fluctuations at zero tempera- 
ture. For simple models there is a correspondence be- 
tween classical and quantum phase transitions such that 
the universal behavior of a _D-dimensional quantum field 
theory corresponds to the critical behavior of a Z) + 1- 
dimensional classical field theory. There are phenomena, 
however, which cannot easily be understood in terms of 
this correspondence, like the nature of the entanglement 
of the ground state wavefunction (see alsc-Si). The en- 
tanglement properties of the quantum wavefunction of a 
device are crucial for determining its suitability as part 
of a quantum computer (see, for instance^). 

Recently, Osterloh et. al.'' discussed entanglement for 
the translationally invariant transversal Ising model in 
one dimension (see also^) . The authors observed that the 
derivative of the concurrence with respect to the coupling 
constant scales according to the Ising universality class 
close to the quantum phase transition. The concurrence 
is a measure of entanglement between only two spin- 1/2 
systemSi^, but also suitable to characterize entanglement 
also between next-nearest neighbors»i The entanglement 
in the transverse Ising and XY models made up of macro- 
scopic (contiguous) subsystems has been discussed ir^i, 
employing the von Neuman entropy as measure of en- 
tanglement. The analysis of the entanglement properties 
near a quantum critical point can be relevant to the anal- 
ysis of many quantum algorithms, as the Hamiltonians 
used to implement them show gapless behavior at some 
point in the computation^. 

We analyze here the entanglement properties of two 
level systems which are either attached to bulk systems 
tuned near a quantum phase transition, or which un- 
dergo a boundary phase transition, as described, for in- 
stance, by the dissipative two level systemiSiii, or the 
Kondo modeli^. These localized phase transitions are 



due to the coupling to a gapless (critical) environmenlii^. 
We will not consider, on the other hand, the entangle- 
ment between qubits as function of their separationiii^, 
even though we distinguish between next and nearest- 
next neighbors. 

We consider in the next section the properties of qubits 
at the boundary of the Ising model in a transverse field, 
i.e., the model studied in'*-^. We analyze how the entan- 
glement between the two last qubits varies as function of 
the distance to the quantum critical point of the model. 
We also allow the values of the couplings at the bound- 
ary to vary. In section III we discuss the entanglement 
properties of two spins attached to a critical reservoir, as 
the coupling between them varies, inducing a boundary 
critical point. We give the main conclusions of our work 
in section IV. 



II. THE ISING MODEL 

We start from the homogeneous Ising model with open 
boundary conditions characterized by the parameter A. 
The two spins at the end are further connected by an 
additional coupling parameter k. The Hamiltonian is 
thus given by 
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where cr^'^ are the x, z-components of the Pauli matrices. 

To solve the modeljiSiii we first convert all the spin 
matrices to spinless fermions. This is done by per- 
forming the well-known Jordan- Wigner transformation 

at = exp ( iTT ^ c]c, j (c, + c|) (2) 
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As usual, we will neglect the term that involves the op- 
erator exp(i7rA/'), N = 'Y^f=ic\ci in order to preserve 
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the bilinearity of the model. An additional Bogoliubov 
transformation then yields (up to a constant) 

N 

T^ = Y1 ^"■^nVn , (4) 
n=l 

N 

Vn = '^{gn,iCi + hn,icl) (5) 
i=l 

where the gn,i, hn^i, and tOn are determined numerically 
(for the general case) . Due to the unitarity of the Bogoli- 
ubov transformation, Eq. (0 is easily inverted to yield 

N 

= ^{9n,i1ln + hn,iT]i) . (6) 
n=l 

A. Concurrence 




FIG. 1: Left hand side: The nearest neighbor concurrence of 
the open boundary Ising model for different locations relative 
to the end as function of A for N — 101. Right hand side: 
The derivative of the concurrence with respect to A. 



We are interested in the reduced density matrix p{i,j) 
represented in the basis of the eigenstates of ■ It is for- 
mally obtained from the ground-state wave function after 
having integrated out all spins but the ones at position 
i and j. As measure of entanglement, we use the con- 
currence between the two spins, C{p{i,j)). It is defined 
as 



C{p{iJ)) = max{0, Ai - A2 - A3 - A4} 



(7) 



where the A^ are the (positive) square roots of the eigen- 
values oi R = pp in descending order. The spin flipped 
density matrix is defined as p = <Ty® ayp*ay (8) cry, where 
the complex conjugate p* is again taken in the basis of 
eigenstates of . It will be instructive to also consider 
the "generalized concurrence" 



C*(p(z,j))-Ai-A2-A3-A4 



(8) 



The reduced density matrix p{i,j) —^p- from now on 
we drop the indices i and j - can be related to correla- 
tion functions. For this, we write the ground-state wave 
function as the superposition of the four states 

m = I m^n) + 1 nmi) + 1 imn) + 1 ii>i<^ii>, 

where the first ket denotes the z-projection of the two 
spins at position i and j. The matrix element P||,xx — 
{4>n\(t>ll)' e.g., is thus given by {'^t<^p, where 

a± = (af ±af)/2.^ 

Due to the invariance of the Hamiltonian under af = 
—erf, at least eight components of the reduced density 
matrix are zero (for finite N). The diagonal entries read: 
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The non-zero off-diagonal entries are 
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P- 
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PUAI 



(13) 
(14) 



The positive square roots of the eigenvalues of R are then 
given by | ^/pipl ± p+ \ and | ^/p2p3 ± P- 1 • Due to the 
semi-definiteness of the density matrix p, we can drop 
the absolute values, i.e., y^piPi ± p+ and ^p2Pz ± P-- 

We now define h = pipA-p2P3 = 4((crf cr|) - (erf) (cr|)) 
and I2 ^ p\ — p^ = —{crfaj) {afa^). For a homogeneous 
model, we have /i > and /2 > Ofii The largest eigen- 
value of Eq. {T)) is thus given by Ai = y/pifH + \p+\ and 
the concurrence reads 



C*(z,j) = 2(|p+|-V^) 



(15) 



We note that the above expression also holds for the gen- 
eralized boundary conditions. For a homogeneous sys- 
tem, it can be further simplified to 

C*(*,j) = (a:,,-l)/2 (16) 

where we introduced the total order Oi j = 

B. Numerical Results 

1. Open boundary conditions 

We first consider the nearest neighbor concurrence of 
the Ising chain with open boundaries (k = 0) for a fixed 
number of sites N = 101 as parameter of A, but for vari- 
ous positions relative to the end of the chain. The results 
are displayed on the left hand side of Fig. ^ As expected, 
the concurrence of the periodic model is approached as 
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FIG. 2: The next-nearest neighbor concurrence of the open 
boundary Ising model for different locations relative to the 
end of the chain as function of A for A'^ — 101. 




FIG. 3: Left hand side: The (generalized) nearest neigh- 
bor concurrence of the closed Ising chain for various coupling 
strengths k as function of A. Left hand side: C*(l,2). Right 
hand side: C(2,3). 



one moves inside the chain and the difference between 
C(50, 51) and C{i, z -I- 1) of the periodic system is hardly 
seen. Nevertheless, the derivative of the concurrence with 
respect to the coupling parameter A, C = dC/dX, still 
shows appreciable differences for A « 1 (right hand side 
of Fig. [TJ. 

We also investigated the scaling behavior of the min- 
imum of C'(l,2), Xmin, for different systems sizes up to 
N — 231. We did not find finite-size scaling behavior for 
the position of the minimum as is the case for the trans- 
lationally invariant model^. The curve of C'(l,2), shown 
on the right hand side of Fig. ^ is thus already close 
to the curve for N —^ oo with a broad minimum around 

Xjnin ~ 1.1. 

The absence of finite-size scaling of the concurrence is 
also manifested in the case of the next-nearest neighbor 
concurrence for different system sizes N. Whereas for 
the periodic system the maximum oi C{i, i + 2) decreases 
monotonically for N —^ oo* there is practically no change 
of C(l, 3) of the open chain for A^>51. 

In Fig. 121 the generalized next-nearest neighbor con- 
currence C*{i,i + 2) of the open boundary Ising model 
is shown for different locations relative to the end of the 
chain as function of A for TV = 101. On the left hand 
side of Fig. |21 results are shown for sites close to the end 
of the chain. Notice that the generalized concurrence 
becomes negative for i = 2, 3 for A > 1 which is not 
related to the quantum phase transition. The crossover 
of the boundary behavior to the bulk behavior is thus 
discontinuous. On the right hand side of Fig. the 
next-nearest neighbor concurrence approaches the result 
of the system with periodic boundary conditions as one 
moves inside the chain. 
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FIG. 4: The generalized next-nearest neighbor concurrence 
of the closed Ising chain for various coupling strengths k as 
function of A. Left hand side: C*(l,3). Right hand side: 
C*(2,4). 



2. Generalized boundary conditions 

We now discuss the concurrence for the generalized 
boundary conditions, introducing the parameter k. On 
the left hand side of Fig. O the generalized concurrence 
of the first two spins C*(l,2) is shown as function of A 
for various coupling strengths k — 0, .., 20A and N — 101. 
For increasing k > 0, the curves indicate stronger non- 
analyticity at A w 1. For k>20X, the generalized con- 
currence becomes negative around A = 1 and is "signif- 
icantly" positive only in the quantum limit of a strong 
transverse field (A < 1). A similar behavior of the con- 
currence is also found in the case of finite temperatures»S, 

On the right hand side of Fig. O the concurrence of 
the second two spins C(2, 3) is shown. All curves display 
similar behavior. There is thus a rapid crossover from 
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the boundary to the bulk-regime and the concurrence of 
periodic boundary conditions is approached for all k as 
one moves further inside the chain. 

To close, we discuss the next-nearest neighbor concur- 
rence C(«,? -I- 2) for various values of n and N — 101. 
On the left hand side of Fig. 01 the generalized concur- 
rence of the first and the third spin, C*(l,3), is shown. 
For K<1, C*(l,3) is positive for all A. For k>1, C*(1,3) 
first becomes negative for A < 1. For k^1.5 C*(1,3) is 
negative for all A. On the right hand side of Fig. ^ the 
generalized concurrence of the second and the forth spin, 
C*(2,4), is shown. For k<A/2, the C*(2,4) is negative 
for A>1. For k>3A/2, the C*(2,4) is negative for A<1. 
Nevertheless, the maximum value is close to A = 1 for all 
cases. 

We finally note that the third neighbor concurrence 
remains zero for all A and all k. 



III. BOUNDARY PHASE TRANSITION 
A. The model. 



cutoff, and a = 1. The Kondo model can be mapped 
onto this modelSS by taking A oc J_l and 1 — a cx . 

To understand the equivalence between these two mod- 
els, it is best to to consider the limit // J ^ 1 (the transi- 
tion takes place fr all values of this ratio). Let us suppose 
that / > so that the Ising coupling is antiferromag- 
netic. The Hilbert space of the two impurities has four 
states. The combinations | IT) and | ||) are almost de- 
coupled from the low energy states, and the transition 
can be analyzed by considering only the | tJ,) and | |t) 
combinations. Thus, we obtain an effective two state sys- 
tem. The transition is driven by the spin flip processes 
described by the Kondo terms. These processes involve 
two simultaneous spin flips in the two reservoirs. Hence, 
the operator which induces these spin flips leads to the 
correspondence A «-> J^/(/ajc). The scaling dimension 
of this term, in the Renormalization Group sense, is re- 
duced with respect to the ordinary Kondo Hamiltonian, 
as two electron-hole pairs must be created. This implies 
the equivalence 2 — a ^ J^. Hence, the transition, which 
for the ordinary Kondo system takes place when changing 
the sign of Jz now requires a flnite value of Jz- 



In order to observe critical behavior of the concurrence 
at the boundary, we thus have to consider a different 
model, i.e., we have to introduce an isotropic coupling 
from spin N to spin 1. This will introduce an interac- 
tion term which contains four fermionic operators and a 
simple solution is thus not possible anymore. 

The model with isotropic coupling is similar to the 
model considered by Garst et. ai, who discuss two Ising- 
coupled Kondo impurities^* (see alscM,). We will consider 
the model studied ini^. It describes to Kondo impurities, 
attached to two different electronic reservoirs, which in- 
teract among themselves through an Ising term. We can 
write the Hamiltonian as: 



H 



— 'HkI + + ISzlSz2 



k 



k.k' -i^LM 



(17) 



We consider the entanglement of the two spins, by writing 
the 4x4 reduced density matrix in terms of the values 
of 5*^1 and Sz2- 

The system described by Ea. ((T7|) undergoes a 
Kosterlitz-Thouless transition between a phase with a 
doubly degenerate ground state and a phase with a non 
degenerate ground state. This transition is equivalent to 
that in the dissipative two level systemiSiii as function of 
the strength of the dissipation. We define the dissipative 
two level system as: 

Utls = ^(t,+Y, \k\b[bk + Xaz ^(^I + bk) (18) 

k k 

The strength of the dissipation can be characterized by 
a dimensionless parameter, a cx A^, and the model un- 
dergoes a transition for A = S/lOc <^ 1, where ujc is the 



B. Calculation of the concurrence. 

The 4x4 reduced density matrix can be decomposed 
into a 2 X 2 box involving the states | ti) and | It), which 
contains the matrix elements which are affected by the 
transition, and the remaining elements involving | H) 
and I II) which are small, and are not modified signifi- 
cantly by the transition. Neglecting these couplings, we 
find that two of the four eigenvalues of the density matrix 
are zero. The other two are determined by the matrix: 



(19) 



where the operator a is defined using the standard no- 
tation of the dissipative two level system, Ea. (|18|l . The 
entanglement can be written as: 



(20) 



The value of (az) is the order parameter of the transition. 
The value of {(Jx), at zero temperature, can be calculated 
from: 

{^.) = ^ (21) 

where E is the energy of the ground state. Using renor- 
malization group arguments, it can be written as: 



E = I 




< a < i 

1 < a < 1 



a - 1 
a > 1 



(22) 
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where C, C and C" are numerical constants. 

If the density matrix is calculated in the absence of 
a symmetry breaking field, {uz) = even in the or- 
dered phase. Then, from Ea. H2UI) . the concurrence is 
given by C = Kfo,)!, which is completely determined us- 
ing Eqs.lini) and In the hmit A/wc < 1 the in- 
teraction with the environment strongly suppresses the 
entanglement. We expect unusual behavior of the con- 
currence for a — \/2 and a = 1. The point a — \/2 
marks the loss of coherent oscillations between the two 
statea^iiSS, although the ground state remains non degen- 
erate. Following the analysis in- , we analyze the behavior 
of dC /da, as a is the parameter which determines the po- 
sition of the critical point. The strongest divergence of 
this quantity occurs for a = 1/2, where: 

On the other hand, near a — 1 the value of dC /da is con- 
tinuous, as the influence of the critical point has a func- 
tional dependence, when a ac, of the type e^'^/'-"^^"^. 
This is the standard behavior at a Kosterlitz-Thouless 
phase transition. This result suggest that the entangle- 
ment is more closely related to the presence of coherence 
between the two qubits than with the phase transition. 
The transition takes place well after the coherent oscil- 
lations between the | ||) and | J,]") states are completely 
suppressed. 

IV. SUMMARY 

We have first calculated the entanglement between 
qubits at the boundary of a spin chain, whose param- 



eters are tuned to be near a quantum critical point. The 
calculations show a behavior which differs significantly 
from the that inside the bulk of the chain. Although the 
spins are part of the critical chain, we find no signs of 
the scaling behavior which can be found in the bulk. We 
use the same approach as done previously for the bulk^, 
although it should be noted that the existence of a finite 
order parameter in the ordered phase will change these 
results if the calculations were performed in the presence 
of an infinitesimal applied field. 

We have also considered the entanglement between two 
spins coupled to a dissipative environment and which 
undergo a local quantum phase transition. The system 
which we have studied belongs to the generic class of sys- 
tems with a Kosterlitz-Thouless transition at zero tem- 
perature, like the Kondo model or the dissipative two 
level system. The most remarkable feature of our results 
is that the entanglement properties show a pronounced 
change at the parameter values where the coherent quan- 
tum oscillations between the qubits are lost, and not at 
the location of the proper phase transition, where the 
ground state becomes degenerate. At this point, how- 
ever, the interaction with the environment has rendered 
the dynamics of the qubits extremely incoherent. 
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